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FOUNDATIONS AND CRITICAL POINT 
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Abstract. The signed loop method is a beautiful way to rigorously 
study the two-dimensional Ising model with no external field. In this 
paper, we explore the foundations of the method, including all details 
that have so far been neglected or overlooked in the literature. We also 
demonstrate how the method can be applied to prove some classical 
results about the Ising model on the square lattice. This in turn leads to 
new explicit formulae for the free energy density and two-point functions 
in terms of sums over loops, valid all the way up to the self-dual point. 
It follows that the self-dual point is critical both for the behaviour of 
the free energy density, and for the decay of the two-point functions. 



1. Introduction and main results 

The Ising model [IlJ is probably the most studied model from statistical 
mechanics, and has been approached using many different ideas and from 
various angles. Most famous are the original algebraic methods of Onsager 
and Kaufman [Hl[T5l[T8], use< ^ by them to compute the free energy and 
study correlations (see also [H]). In the fifties, Kac and Ward [12] proposed 
a very interesting alternative method for approaching the Ising model, based 
on studying configurations of signed loops. This method is combinatorial 
in nature, and as such often referred to as the combinatorial method. How- 
ever, in order to avoid confusion with other approaches which also include 
combinatorial aspects (for instance, Kasteleyns dimer approach as exposed 
in [17]), we prefer to refer to the method as the signed loop method. 

Over the years, a number of articles developing the Kac-Ward method 
have appeared in the physics literature, of which the most relevant are [31E1 
[2T j l22 t [25]. However, from a mathematical point of view, these papers, as well 
as the recent review [5], leave a lot to be desired in terms of mathematical 
rigour and technical details. Moreover, doubts have been cast on the very 
validity of the whole method to begin with, not only in the years following 
the Kac-Ward paper, but still relatively recently by Dolbilin et al. jS], who 
rightly pointed out an error in Vdovichenko's article [25] (reproduced in [161 
§151]) on the signed loop method. 
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With the present paper, we aim to remove these doubts and deficiencies 
once and for all, by providing complete, rigorous and detailed proofs of the 
combinatorial identities central to the signed loop method, in great gener- 
ality. In particular, we show that the error in Vdovichenko's article can be 
corrected by including an additional factor (the multiplicity defined below) 
into the loop weights. 

The main objects of study in the signed loop method are graph generating 
functions. As we shall see, the one-point and two-point functions of the Ising 
model on 1?, as well as the partition function, can all be expressed in terms 
of graph generating functions. As a result, the signed loop method can 
be used to prove both classical and new results for the Ising model. In 
particular, we show that the method leads to a proof of the fact that the 
self-dual point is in fact critical, both in terms of spontaneous magnetization, 
and as the only point where the free energy density is not analytic. 

Needless to say, there are other approaches to the Ising model which 
appeal to mathematicians, and to probabilists in particular. We mention 
the approach of Aizenman, Barsky and Fernandez pQ, who prove sharpness 
of the phase transition using differential inequalities. We also mention the 
fermionic observables originally introduced by Smirnov [24J to study the 
Ising model at criticality, but which have been used in [2] in an interesting 
way for yet another derivation of the value of the critical point. 

This paper is organised as follows. We present our key identities con- 
cerning graph generating functions and signed loops in Section and we 
present our results for the Ising model on 1? in Section 11.21 The proofs of 
our combinatorial identities are given in Section [21 and the proofs of our 
results for the Ising model can be found in Section [3l 

1.1. The combinatorial identities. Although our applications are in Z 2 , 
we will occasionally have to augment the graph on which we work. Therefore, 
it will be necessary to study the signed loop method on a more general graph 
structure, which is also of interest in itself for further applications. Hence, 
our starting point is a finite graph G = (V, E), embedded in the plane, with 
vertex set V and edge set E. We identify G with its embedding. We assume 
G does not have multiple edges, but we do not assume that G is planar. 
For convenience, although this is not strictly necessary, we assume that the 
embedding is such that edges are straight line segments, and that except for 
the vertices at the two endpoints, no other vertices lie on an edge. 

For practical reasons, we will think of edges as unordered pairs of vertices 
{u,v}, u,v £ V. For brevity, we will also write uv instead of {u,v} for the 
edge between u and v (note that uv = vu). Since we do not assume that G 
is planar, it is possible that two edges intersect in a point which is not a 
vertex. In this case, we say that the two edges cross each other. 

We call a subset F of E even if every vertex in the subgraph (V, F) of G 
has even degree (the empty set is also even). Our main object of study in 
this section will be a type of generating function for these even subgraphs. 
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To define this generating function, we assume that for every edge uv E E, 
we are given an edge weight x uv G R. We write x = {x uv ) uv( zp for the 
vector of all edge weights. Furthermore, we denote by Cp the total number 
of unordered pairs of edges in F that cross each other. The generating 
function Z(x) of even subgraphs of G is then defined as 

(1.1) z( X ) = £ (-ip n x ™- 

even FcE uv&F 

If the graph G is planar, then we can embed G in such a way that there 
are no edges that cross each other, so that Cp = for all F. Note, however, 
that Z(x) depends on the embedding, and that two different embeddings 
of the same (abstract) graph can lead to different functions Z(x). Since we 
identify G with its embedding, this last fact does not concern us here. 

The signed loop method is a way to compute the generating function Z(x) 
by expressing it in terms of a sum over all possible loops in the graph G. 
Loosely speaking, a loop t is a non-backtracking path passing through ver- 
tices Vo, vi, . . . , v n = Vq. It has a sign sgn(^) which is —1 if its total winding 
angle is an even multiple of 2tt, and +1 otherwise, and a multiplicity m(£) 
counting how often the loop "repeats" itself (precise definitions will be given 
below). We associate with each loop i = («o, • • • ,v n -i) a signed weight 

(1.2) ^)=f^n^ +1 > 

and we claim that with this choice of loop weights, the following holds: 

Theorem 1.1. For uv £ E, let d uv denote the maximum of the degrees of 
u and v in the graph G. If \x uv \ < [d uv — 1) _1 for all uv £ E, then 

(1.3) Z(x) = exp(j2 w ( i > x ))' 

H inG ' 

Under the condition of Theorem 1 the loop weights are absolutely sum- 
mable, but we shall see that under a significantly weaker condition, we still 
obtain an equality similar to (II. 3D if we sum the loop weights in order of 
increasing length of the loops. Before we explain this, we will first give the 
precise definitions of loops in G, and of their signs and multiplicities. 

We start by introducing the notion of a path, since a loop is a specific 
type of path. A path of n steps in G is a sequence (vo,v±, . . . ,v n —i) £ V n 
such that ViVi + \ € E for every i satisfying < i < n — 2, and v «+2 / Vi for 
every i such that < i < n — 3 (paths are non-backtracking). If the path 
(vo, . . . ,v n -\) satisfies vo^n-i £ E, vq ^ « n -2 and v\ / v n -\, then it is also 
called a rooted loop of n steps, rooted at vo. In this case, all rotations of the 
sequence (vo, . . . ,f n -i) are also rooted loops of n steps, and so are all the 
rotations of the reverse sequence (v n -\, v n -2, ■ ■ ■ , vq). 

We think of all these rooted loops as alternative representations of the 
same unrooted loop, and we want to choose one of these representations 
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Figure 1. The turning angle from the vector v — u to the 
vector w — v (left). The loop (v%, vi, ^3, V4,v$,vq) on the right 
has sign —1, the loop (vi, V2, V3, v$, V4, Vq) has sign +1. 

as the one identifying the unrooted loop uniquely. To do so, we order the 
vertices of G lexicographically by their coordinates in the plane. Then there 
is a unique sequence (vq, . . . , t> n -i) which is lexicographically smallest among 
the collection consisting of all its rotations, and all rotations of the reverse 
sequence (i> n — 1, • • • ,vo). Henceforth, when we speak of a loop (without the 
qualification rooted), we shall mean such a sequence which is smallest among 
the corresponding collection of permutations. 

If £ = (vq, . . . , u n _ 1) is a rooted or unrooted loop, we shall make the 
identification Vj = Vj mo d „ for all j 6 Z. We say that a loop t is edge- 
disjoint if ViVi + \ 7^ VjVj + \ for all i,j £ {0, ...,n — 1} with i 7^ j. If £ 
is not edge-disjoint, it might be the case that the sequence (vq, ... 
is periodic, in which case we call £ a periodic loop. The multiplicity of £, 
denoted by m(£), is its number of steps divided by its smallest period. In 
particular, the multiplicity of every nonperiodic loop is 1. 

The sign of a loop is defined in terms of its winding angle, as follows. 
Given two distinct edges uv and vw, we define Z(v — u, w — v) as the turning 
angle in the plane from the vector v — u to w — v, see Figure Q] (left). The 
winding angle a(£) of a loop £ = (vq, . . . ,u n -i) is simply the sum of all 
turning angles along the loop, that is, 

n— 1 

(1-4) a{£) = Z(u i+ i - vi, v i+2 - v i+1 ). 

We now define the sign sgn(£) of £ as 

(1.5) sgn(*) = -expQa(*)). 

Observe that the winding angle of every loop is a multiple of 2tt (here we 
use the fact that the edges of G are straight line segments), hence the sign 
of a loop is either +1 or —1. 

Remark. If a loop is edge-disjoint, it follows from Whitney's formula [26] 
that the sign of the loop is —1 if the loop crosses itself an odd number of 
times, and +1 if it crosses itself an even number of times (see Figured]). For 
loops that are not edge-disjoint, it may not be so clear what is meant by a 
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Figure 2. A square lattice wrapped around a torus (right) 
and a representation of it in the plane (left). The gray square 
corresponds to the torus, the dotted lines and open circles are 
the additional edges and vertices used for the embedding. 



"crossing", but definition (|1.5|) makes sense for both kinds of loop. Readers 
familiar with the recent work of Smirnov will notice the similarity between 
definition (|1.5|) and the definition of the fermionic observable in [21], which 
is also in terms of winding angles. 

Note that there are infinitely many loops in G (unless G has no cycles), 
since loops can have any number of steps. In what is to follow, it becomes 
essential to sum over loops in the correct order; we will collect loops of a 
given length together, and sum in order of increasing length. However, we 
will not always identify the length of a loop with its number of steps. 

For instance, in our applications to the Ising model on Z 2 , we will have 
reasons to consider loops on a square lattice wrapped on a torus. However, 
for the signed loop method and our Theorems 11.11 and 11.21 it is essential 
to work on a graph G embedded in the plane. Instead of working on the 
torus directly, we must therefore work on a suitable representation of it in the 
plane. To keep the computations simple, the most convenient representation 
will be a rectangular piece of 1?, to which some extra edges and vertices are 
added to connect opposite sides, as illustrated in Figure [21 

Some of the edges in the constructed graph G now do not correspond 
to edges that can be traversed by a loop on the torus, and therefore, steps 
taken along these added edges should not contribute to the length of a loop. 
We call these extra edges additional edges, and denote the set of additional 
edges by Ea- The edges in E \ Ea are called representative edges, since 
they do correspond to edges in the represented graph. In general, these 
considerations lead us to allow the edge set E of the graph G we work on to 
be divided into a set Ea of additional edges and a set E\Ea of representative 
edges. For reasons that will become clear, we must assume that the set Ea 
is such that the graph (V, Ea) is free of cycles. 

Additional edges will also play another role. In our study of correlations 
in the Ising model, they are used to artificially connect two spin variables 
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whose correlation we want to study. Again, also in this role, we do not want 
the additional edges to contribute to the length of a loop. 

For these reasons, we will measure the length of a loop in terms of the 
number of representative edges it traverses, ignoring any steps taken along 
additional edges. To be precise, we formally define the length r{i) of the 
loop I = (vq, . . . , v n -i) as the number of i in {0, . . . , n— 1} such that UfUi+i G 
E\Ea- Note the distinction between the length of a loop and its number of 
steps. By C r we denote the collection of all loops of length r, and we write 
C = U^Li &r for the collection of all loops. We can now write 

(1.6) f r {x)=Y J w(i;x) 

for the sum of weights of all loops of length r. Under the condition of 
Theorem II. 1\ the loop weights are absolutely summable. This implies in 
particular that Z{x) equals exp(V^ r f r (x)), but we shall see that the latter 
policy of summing the loop weights by length can already be used under a 
significantly weaker condition. 

To formulate this condition, we need to introduce the transition ma- 
trix A(x). If uv is an edge of G, then by ui we will denote the directed 
edge from u to v. The matrix A(x) will be indexed by the directed represen- 
tative edges of G. Given two directed representative edges ut> and wt, we 
say that v is linked to w if either v = w, or there exists a sequence of distinct 
additional edges «ii>2, ^2^3, . . . , v n -\v n such that v = v\ and v n = w. In the 
former case, if v = w and u ^ z, we write 

Z(u1), wz) = Z(v — n, z — w) 

for the turning angle from ui to wt. In the latter case, the sequence 
{v\, . . . , v n ) is a path (the chain) linking v to w, passing through additional 
edges only, and we say that u v w via (v±, . . . ,v n )". By our assumption 
that the additional edges form no cycles, there can be at most one such path. 
Hence, without ambiguity, if v is linked to w in this way, we can define 

Z(ui, wt) = Z(f — U,V2 — v) 
n-2 

+ Z(v i+ i - Vi,V i+2 - V i+ i) + Z(w - V n -!, Z - w). 
i=l 

The transition matrix A(x) is defined as follows. Write A^^(x) for the 
entry of the matrix with row index ui and column index wt. Then 

if v = w and u 7^ z; 
if v w via (y±, . . . , v n ); 
otherwise. 



(1-7) A^^(x) = ^ 



XitfiC 



iZ(ui,wi)/2 



n—1 , v 

rp FT , e iAuv,wl)/2 

i=l 





SIGNED LOOP APPROACH TO THE ISING MODEL 



7 



Let \i(x), i = 1,2,..., 2\E \ Ea\, denote the eigenvalues of A(x). We will 
show that if maxj|Ai(x)| < 1, then the f r (x) are absolutely summable. This 
leads to an identity which forms the core of the signed loop method: 

Theorem 1.2. If maxj|Aj(a;)| < 1, then 

(1.8) Z(x) = exJjrfrix)) = ^det(l-A(:r)). 

\ r=1 / 

An alternative proof of the (Kac-Ward) formula Z(x) 2 = det(l— A(:r)) for 
the case of a planar graph G with straight edges, and without consideration 
of the possibility of including additional edges, was given in [7]. Our results 
are not only more general, but also express Z{x) in terms of loops and 
loop configurations. Both these facts are used in an essential way in our 
applications to the Ising model. Extensions of the Kac-Ward formula to 
graphs embedded in surfaces of higher genus are studied in [I] . 

1.2. Applications to the Ising model. In this section, we consider ap- 
plications of Theorem ll.2l to the two-dimensional Ising model. We focus our 
attention here on studying the phase transition of the isotropic Ising model 
on the square lattice 1?, in terms of the behaviour of the free energy density 
and of the two-point functions. Similar results as the ones we obtain here 
for the square lattice can also be obtained for the triangular and hexagonal 
lattices with the signed loop method. In fact, the method applies to the 
Ising model on even more general graphs, and allows one to study general 
fc-point functions. We will go into these issues in a subsequent paper. 

In our applications to the Ising model, ||-|| will denote the Li-norm in T? 
or R 2 . Thus, ||u — i>|| measures the graph distance between vertices u and v 
of 1? or Z 2 *, the dual of 1?. The notation ||-|| is used for the maximum 
norm on the relevant space, and o = (0, 0) denotes the origin in 1? . 

We shall start with a brief review of the definition of the Ising model and 
of the low- and high-temperature expansions, through which the connection 
with the graph generating function Z(x) of Section 11.11 will be made. We 
consider the model on finite rectangles in the square lattice 1?. We call a 
graph G = (V, E) a rectangle in 1? if the vertex set V consists of all points 
of Z 2 contained in a rectangle [a, b] x [c, d] in R 2 , and for u, v S V, the edge 
set E contains the edge uv if and only if ||u — v\\ = 1. We say that vertices 
u and v are neighbours if uv G E, and we define the boundary dG of G as 
the set of vertices that have less than 4 neighbours in G. By G* = (V*,E*) 
we shall denote the weak dual graph of G, i.e. the rectangle in 1?* whose 
vertices are the centres of the faces of G (see Figure El left). 

To define the Ising model on G, we associate with each vertex v G V a 
spin variable o~ v , taking values in {— 1,+1}. Thus, the space of all possible 
spin configurations is Q = {— 1,+1} V . Sometimes we shall impose positive 
boundary conditions, by which we mean that we restrict ourselves to the set 
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of spin configurations 

Q.+ = {a € n : o u = +1 if u £ dG}. 

In contrast, when we speak of free boundary conditions, we work with the 
unrestricted set of spin configurations fi free = Q. 

To each spin configuration a = (cr v ) ve v we associate an energy given by 
the Hamiltonian 

H(a) = - a u a v . 

uv£E 

If we interpret the spin variables as magnetic moments, this Hamiltonian 
describes ferromagnetic interactions, since the energy is lower the more pairs 
of neighbouring spins are in the same state. According to the formalism of 
statistical mechanics, the equilibrium probability distribution of the spin 
configurations at the inverse temperature (3 is now given by the Gibbs- 
Boltzmann distribution 

(1.9) If {a) = ±e-f>*V = ± II ^\ a E , 

P P uv&E 

where □ £ {free, +} stands for the imposed boundary condition, and Zg, 
the partition function, is given by 



(i.io) zf = y: = e n < 



1.2.1. Low- and high-temperature expansions. The partition function of the 
Ising model is closely related to the graph generating function Z{x) from 
Section 11.11 To see this, we now consider the so-called low-temperature 
and high-temperature expansions. More details on these expansions and the 
related duality of the Ising model can be found in [23, Section II. 7]. For our 
purposes, the low-temperature expansion is best considered in the case of 
positive boundary conditions. It is not difficult to see that in this case, there 
is a 1-1 correspondence between the even subgraphs of the weak dual G* and 
the spin configurations in Q + : given a £ one obtains the corresponding 
even subset F(a) of E* by including the edge dual to uv in F(a) if and only 
if a u 7^ a v , for every uv € E. See Figure [3] (left) for an illustration. 

Note that by this correspondence, if F C E* is even, then every edge in F 
separates two spins that have opposite sign. Since the energy increases by 2 
if cr u a v changes from +1 to —1, this means that every edge uv € F comes 
with a "cost" of exp(— 2/3). It follows that we can write 

p; w = «!) jj „-» 



where 

(1.12) Z+ = ezptfm £ J] 

even FcE* uv£F 



e~^. 
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Figure 3. Left: the graph G and its weak dual G*, with an 
even subgraph of G* marked by bold dashed edges. The spins 
in the gray squares have value —1, the rest have value +1. 
Right: the subgraph of G drawn with bold edges contributes 
cr x cr y a z 'u a v (tanh/3) 5 in the high-temperature expansion. 



This is the low-temperature expansion of the partition function for posi- 
tive boundary conditions. Observe that up to the factor exp((3\E\), this ex- 
pansion takes exactly the form (jl.ip of the graph generating function Z(x) 
for the dual graph G* (in which no edges cross each other), if we set the 
edge weights of all dual edges uv 6 E* equal to x uv = exp(— 2/3). 

We now turn to the high-temperature expansion, for which we impose 
free boundary conditions. The expansion will be over even subgraphs of 
the graph G, rather than of the dual G*. Unlike in the low-temperature 
expansion, these subgraphs do not have a clear geometric interpretation. 
We will therefore take some time to explain how they arise. 

The high-temperature expansion starts from equation (jl.lOp and the ob- 
servation that a u a v can only take the values —1 or +1. Since exp(±/3) = 
cosh f3 ± sinh /3, we see that 

4 ree = (cosh /3)l s l ^ II i 1 + a * a « tanh /3 ) • 
cren frec uvge 

The next step is to expand the product over uv G E. Each term in the 
expansion will be a product of factors obtained by choosing for each edge uv 
whether 1 is taken as a factor, or a u a v tanh j3. The expansion can therefore 
be written as the sum of the products of factors over all choices one can 
make for the factors coming from each edge uv. We can represent each 
choice graphically by removing the edge uv if we choose the factor 1 for this 
edge, and keeping uv if we choose the factor a u a v tanh f3. This gives a 1-1 
correspondence between all terms in the expansion, and all F C E (not just 
the even ones). See also Figure [3] (right). 

Using this correspondence, and then interchanging the order of summa- 
tion over a and F, we may now write the partition function as 

^(cosh/?)^ £ n^ deg(u,F) 
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where deg(u, F) denotes the degree of u in the graph (V, F). Note that the 
sum over a vanishes unless deg(u, F) is even for all u E V, in which case the 
sum yields simply 2' v \. Therefore, 

(1.13) 4 ree = 2l y l(cosh/?)l £; l I] tanh ^' 

even FcE uv£F 

Again, up to a multiplicative constant, the expansion takes exactly the 
form (II. ip of the graph generating function Z(x), this time for the graph G, 
if we set the edge weights equal to x uv = tanh/3. 

1.2.2. Results for the Ising model. The connection between the low- and 
high-temperature expansions and the generating function Z(x) can be ex- 
ploited to derive a number of results for the Ising model. In this paper, we 
focus on specific aspects of the phase transition in the isotropic Ising model 
on the square lattice. We shall study in particular the free energy density 
and the behaviour of the two-point functions. 

In order to apply Theorem II. 2^ we must obtain a bound on the spectral 
radius maxj|Aj(x)| of A(x). Our following result therefore provides the key 
to applying the signed loop method to the Ising model on 1?: 

Theorem 1.3. For a finite rectangle G in 1? with no additional edges, 

max|Ai(x)| < (y/2 + ljIML, 

i 

and |/ r (x)|<2|F|(v^+l) r |NI«. 

Recall that we will take our edge weights to be respectively exp(— 2/3) 
and tanh/3 in the low- and high-temperature expansions. Thus, Theorem ll.3l 
says that we can apply the signed loop method for all temperatures up to 
the self-dual point, which is given by the value of /3 C satisfying 

(1.14) exp(-2/3 c ) = tanh/3 c = V2 - 1. 

We write /3 C here for the self-dual point, because our results actually imply 
that /3 C is the critical point of the Ising model in two senses: it is the only 
point where the free energy density is not analytic, and also the point where 
the behaviour of the two-point functions changes from staying bounded away 
from zero, to decaying exponentially with distance. A computation of the 
latter interpretation of the critical point was recently also performed in [2], 
using a very different method based on Smirnov's fermionic observable [24] . 

The reader may be surprised to learn that the signed loop method applies 
all the way up to the critical point, while classical arguments like the Peierls 
and Fisher arguments [5J 10. 2(1] stay far from it. What happens here is the 
following. If we simply count the number of non-backtracking (unsigned) 
loops in 1? of n steps starting from a vertex v, this number grows like 3 n . 
However, if we take the signs of these loops into account, first counting all 
loops that have positive sign, and then subtracting those that have a negative 
sign, the obtained number only grows (in absolute value) like (V2+l) n . This 
reduction in growth rate is precisely enough to take us to the critical point. 
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We now turn to our results for the Ising model. First we consider the 
Helmholtz free energy F^, which is related to the partition function through 

Ff = ~^- 

The free energy is an important quantity, because it encodes all kinds of 
thermodynamic information about the system: derivatives of the free en- 
ergy with respect to temperature are related to physical properties such as 
internal energy and specific heat. In particular, the formalism of statistical 
mechanics predicts that phase transitions coincide with discontinuities or 
other singularities in derivatives of the free energy. However, since the free 
energy of any finite system is always an analytic function of the parameters, 
such singularities can only occur in the infinite- volume limit. 

We will therefore consider limits as our rectangles G = (V, E) tend to 1? 
(i.e. such that every finite subgraph of 1? is eventually contained in G). 
Since now we no longer work on a fixed graph, we shall make the depen- 
dence on the graph explicit in our notation by writing Zq ^ for the partition 
function of the Ising model on G, and Fq a for the free energy. It is a 
well-known fact that the infinite- volume limit of the free energy density, i.e. 

lim -i-F5„= lim — lnZ9 R =: f(B), 

exists (which we will prove). It is also not difficult to show that the limit is 
the same for all boundary conditions, see for example [23} Section II. 3]. 

We will use Theorem 11.21 to derive an explicit power series expression 
for f(/3). This implies directly that f((3) is an analytic function of f3 on 
(0,/3 c ) U (/3 c ,oo), and moreover, it allows us to derive the formula for /(/?) 
which was first obtained by Onsager in [18] : 



Theorem 1.4. The free energy density /(/3) is an analytic function of ft on 
(0, j3 c ) U ((3 c ,oo), where it is given by the Onsager formula 

— — — 7T / / In [4 cosh 2 2/3 — 4sinh2/3(cos oj\ + cosw2)l dcoi du)2- 

The functions / and u = d(f3f)/d(3, which is the internal energy density 
of the system, are both continuous functions of (3. However, Onsager has 
shown in [18J that the derivative of u with respect to temperature, which is 
the specific heat, diverges as /3 approaches j3 c . This shows that fi c is critical 
for the behaviour of the free energy. But we would really like to show that 
/3 C is also critical for the magnetic behaviour of the Ising model. Next, we 
therefore consider the one-point and two-point functions above and below f3 c . 

To simplify the notation, and following the physics literature, we write 

(1.15) (o"«)g,/3 = E GA au ^ ( cr «^)B,/3 = E GA au(7 v)' 

where Eq a denotes expectation with respect to the Ising measure Pq a- Our 
next result says that for positive boundary conditions, at all temperatures 
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below the critical point, the system has a net magnetization per unit volume 
in the infinite- volume limit: 

Theorem 1.5. For all f3 £ (/3 c ,oo) and fixed u, v 6 1? , the limits 

lim {<r v )%p =: (<t„)+ „ and lim (a u a v ) % p =: (a u a v )+ 2 „ 

both exist and are strictly positive. Furthermore, we have that 
lim {a u o- v )+ 2 = [{a )+ 2 J 2 > 0. 

\\u— v\\ — >oo ' K ' r 

For free boundary conditions, by symmetry under changing the values 
of all spins, (<J u }q^ = 0, so the infinite- volume limit {o~ u )^ 2 e p trivially ex- 
ists and is 0. The next theorem addresses the behaviour of the two-point 
functions at high temperatures: 

Theorem 1.6. For all (3 £ (0, j3 c ) and fixed u,v G ^? , the limit 

exists. Moreover, we have that 

„ I™ (a u ^)l^=[(a o )|!y 2 = 0. 

\\u— v\\— >oo 

Our proofs show that in all cases in Theorems 11,51 and 11.61 the conver- 
gence is exponentially fast. Note the contrasting behaviour of the two-point 
functions above and below /3 C : at low temperatures they stay bounded away 
from 0, while for high temperatures they decay to 0. This is an indication of 
the phase transition of the Ising model in terms of its magnetic behaviour. 

However, we have used different boundary conditions above and below the 
critical point. Ideally, we would like to show that there is no magnetization 
in the high-temperature regime even for positive boundary conditions, and 
that at low temperatures, the two-point functions stay bounded away from 
zero for free boundary conditions as well. We believe that this can be shown 
using the signed loop method, but leave these issues for a subsequent paper. 

2. Proofs of the combinatorial identities 

We now turn to the proof of our main Theorems 11.11 and 11.21 Recall that 
here G = (V, E) is a general finite graph embedded in the plane, poten- 
tially containing crossing edges or additional edges. The proof proceeds in 
a number of steps. In the first step, detailed in Section \2.1\ we will identify 
each even subgraph (V, F) of G with a number of edge-disjoint collections of 
loops, and show that the sum of their weights yields precisely the contribu- 
tion of F to the generating function Z(x). In the second step, in Section [2]2j 
we will explain the conditions under which we can express Z(x) in terms of 
Ylr fr( x ) an d det(l — A(x)), under the assumption that the weights of all 
remaining configurations of loops in G of total length r cancel each other. 
The proof of this assumption, the last step, will be carried out in Section l2~3l 
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2.1. Expansion into collections of edge-disjoint loops. We will be con- 
cerned with crossings of loops and paths, and we need to carefully establish 
the relevant definitions first. More specifically, we will consider collections 
{£\, . . . ,£ s } of loops with the properties that all loops £i,. . . ,£ s are edge- 
disjoint, and no two loops in the collection visit a common edge. We call 
these edge- disjoint collections of loops. Intuitively it may be clear what we 
mean by a crossing of such edge-disjoint loops, but some care is needed, so 
we will now give the precise definitions. 

First, consider two paths (u, v, w) and (x, y, z) in G, and let A be the 
union of the two half-lines {v + t(u — v ) : t > 0} and {v + t(w — v) : t > 0}. 
We say that (u,v,w) and (x,y,z) cross each other at the vertex v if v = y 
and the vertices x and z do not lie in the same infinite component of the 
complement of A in the plane. 

Now let l\ = (uq, . . . , u n -i) and £2 = (yo, . . . , u m _i) be two loops that 
form an edge-disjoint pair {£±,£2}- By Cy(£i,£2) we denote the number 
of pairs where < i < n and < j < m, such that the paths 

(tti-i, Ui, Ui + i) and (vj-i,Vj,Vj + \) cross each other at Ui. We call Cy^i,-^) 
the number of vertex crossings between £\ and £2- Similarly, we define 
the number of edge crossings between £\ and £2, denoted Ce(£i,£2), as the 
number of pairs (i,j) such that < i < n and < j < m, and the edges 
UiUi + i and VjVj+i cross each other in G. 

We also need to formally define the number of times a loop crosses itself, 
so consider an edge-disjoint loop £ = (vo, ■ ■ ■ ,v n -i). We define the number 
of vertex self-crossings of £, denoted Cy{£), as the number of pairs 
where < i < j < n, such that Uj, i>i+i) and (yj-i,Vj,Vj+i) cross 

each other at the vertex Vj. The number of edge self-crossings of £, denoted 
Ce(£), is defined as the number of pairs such that < i < j < n, and 
the edges fjfi+i and VjVj+x cross each other in the graph G. 

As we have already mentioned in the introduction, Whitney's formula [26J 
says that the sign of an edge-disjoint loop is —1 if the loop crosses itself an 
odd number of times, and +1 otherwise. In other words, we have 

sgn(£) = (-lfvd)+c E (i) 

if £ is edge-disjoint. We now simply define the sign of an edge-disjoint 
collection of loops {£1, . . . ,£ s } as 

s 

(2.1) sgn{4,...,4} = Hsgn(£ i ) = (-1)£? =1 {<M4)+C^)}_ 

i=l 

If F C E is even, we can decompose F into an edge-disjoint collection of 
loops in such a way, that the union of all edges traversed by the loops is F 
(one way to find such a decomposition is given in the proof of Proposition ^. II 
below). This decomposition is in general not unique. We write T>(F) for 
the set of all possible edge-disjoint decompositions of F, and recall that Cf 
denotes the number of unordered pairs of edges in F that cross each other. 
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Vi 



' v 5 



•v 6 



Figure 4. The neighbours vi,V2, ■ ■ ■ of a vertex v in an even 
subgraph (V,F) are ordered in a clockwise fashion around v. 
Two neighbours that are connected by edges drawn in the 
same line style are paired to each other (see the text). 

Proposition 2.1. For all even subsets F of E we have that 

Y sgn{£ 1 ,...,£ s } = (-lf". 
{t u ...,e s }ev(F) 

Proof. Let {£\ , . . . , £ s } be an edge-disjoint collection of loops. Since l\ , . . . , £ s 
are all closed loops in the plane, any two distinct loops £j and £j from the 
collection necessarily cross each other an even number of times. That is, 
C v (£ h £j) + C E {£ h £j) is even for all i ^ j. Therefore, 

sgn{4, . . . ,4} = (-i)Ei< I < s {^v(4)+c E (4)}+Ei< 1<J < s {Cv(^/ J )+c B (4^)} > 

Furthermore, if {£±, . . . ,£ s } G T>(F), then clearly the total number of edge 
crossings occurring among the loops must coincide with Cp, that is, 

c F = c dti)+ Y c E(ti,tj)- 

l<i<s l<i<j<s 

Hence, it suffices to prove that 

(2.2) Y (-l)£i<i<s c ^)+£i<i<i<3 C ^ £i A) = 1. 

{lu.. .,i s }£V(F) 

Let Vf denote the set of vertices in V whose degree in the subgraph (V, F) 
is nonzero, and let deg(v,F) denote the degree of v in (V,F). For v 6 Vf 
of degree 2k, write v\, . . . , V2k for the endpoints of the edges in F that are 
incident to v. Assume that these vertices are ordered in a clockwise manner 
around v, starting from the lexicographically smallest one (see Figure H|). 
Denote by V v (F) the collection of partitions of the vertices v\, . . . , V2k into 
sets of size 2. We call these partitions pairings at the vertex v. We write 

V(F) = J] V V {F), 

v&V F 



and call an element of V(F) a pairing associated with the subgraph (V,F). 
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We have a natural 1-1 correspondence between V(F) and T>{F). Indeed, 
starting from any vertex v G Vf and any i G {1, . . . , deg(v, F)}, the pairing 
7r G T'(-F) defines a unique rooted loop (uq, ■ ■ ■ ,u n -l) with the properties 
that Uq = Vi, u± = v, and for all j, Uj_i is paired with Uj + \ at the vertex Uj. 
Continuing this way, and replacing each rooted loop obtained by its unrooted 
counterpart, yields an edge-disjoint collection . . . ,£ s } G F>(F). It is easy 
to see that this defines a bijective relation between V(F) and T>(F). 

Using this bijection, we can express the sum in (|2.2p equally well as a sum 
over all pairings. More precisely, for tt G V(F) let ir v denote the pairing it 
induces at the vertex v G Vf, and write C v (ir v ) for the number of crossings 
at the vertex v introduced by this pairing. We call n v even (odd) if C v (tt v ) 
is even (odd). Note that (|2,2p is equivalent to 

ireV(F) v&V F tt v €Vv(F) 

from which we see that it suffices to prove that for all v G Vf, the number 
of even pairings tt v exceeds the number of odd pairings tt v by 1. 

We prove this by induction on the degree 2k of v. Write and for 
the numbers of even, resp. odd, pairings for v of degree 2k. For k = 1 we 
clearly have = 1 and = 0. Now let k > 1, and suppose that we pair 
the vertex v% with Vi at v. Next pair the remaining 2k — 2 neighbours Vj 
of v in all possible ways. For even i, there is an even number of j in between 
1 and i, and therefore the pairing we obtain will be even if and only if the 
pairing of the remaining 2k — 2 vertices is even (see Figured]). Likewise, 
for odd i, the obtained pairing will be even if and only if the pairing of the 
remaining vertices is odd. Since we have k even values for i, and k — 1 odd 
values, this gives N+ = fcJV^. 1 + (fe-l)JV^. 1 and = kN^ + ^k-^N^. 
Hence by the induction hypothesis, — = 1. □ 

From Proposition 12. 1 1 we will now obtain our first main result. Recall that 

z{x) = (-!) cf n 

even FcE uv&F 

The case F = is treated separately: by convention it contributes 1 to the 
sum. Hence, Proposition 12.11 implies that 

Z(x) = l+ J2 sgn{h,...,£ s } Yl x uv . 

even FcE: U 1 ,...l s \eT>(F) uv&F 

Recall that the multiplicity m{€) of an edge-disjoint loop t is 1. Therefore, 
using ()2.ip and the definition f)1.2|) of the weight of a loop, we can write 

s 

z(x)=i+ e n^x). 

even FcE: {h,...,l s \eD(F) i=l 
F+$ 
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Since this is a finite sum, we do not need to worry about the order of 
summation, so we have 

oo s 

Z(x) = i + J2 E E n«(4;a0. 

r=l even FcE: (£i,...,£„}€D(F) i=l 
\F\E A \=r 

If we now denote by T> T the set consisting of all those edge-disjoint collections 
of loops . . . ,£ s } for which the total length Ylt=i r (^i) 1S r > we see that 
we have established the following theorem: 

Theorem 2.2. 

oo s 

z(x)=i+j2 e n^- x )- 

r=l {£i,...,£ s }eV r i=l 

2.2. Extension to all loop configurations. In Theorem 12.21 we have 
expressed the generating function Z(x) as a sum over all edge-disjoint col- 
lections of loops in G. In this section, we will see that if the edge weights are 
sufficiently small, we can drop the condition that the loops have to be edge- 
disjoint, and sum instead over all possible loop configurations in G. Here, a 
loop configuration is simply an ordered sequence {£\, . . . ,£ s ) of loops; there 
is no condition that loops have to be edge-disjoint, nor that two loops in 
the configuration have to be distinct (i.e. it is allowed that £j = £j for some 
i j, which is why we work with ordered sequences of loops now). 

Write C r for the collection of all loop configurations (£i, . . . ,£ s ) satisfying 
r{£\) + ■ ■ ■ + r(£ s ) = r. Some of these loop configurations will consist of 
distinct loops that together form an edge-disjoint collection of loops. Let 
C* denote the subset of C r containing only these edge-disjoint loop configu- 
rations. Observe that if {£\, . . . ,£ s } is an edge-disjoint collection of loops, 
then the corresponding loop configuration (£i, . . . ,£ s ) has s\ permutations. 
Therefore, by Theorem 12.21 we already have that 

CO oo I s 

s(*) = i + EE E 

r=ls=l(£ 1 ,.,{,)eC; "i=l 

but we claim that here we may sum over C r instead of C*: 
Theorem 2.3. 

OO CO 1 S 

z(x) = i + EE E ^tIN*;*)- 

r=l s=l {£ ly ..,£ a )eCr i=l 

Clearly, since C r is a finite set for every fixed r, this result is an immediate 
consequence of the following proposition: 

Proposition 2.4. For all r > 0, 

CO I s 

e e iT n^^)=°- 

s=l (£ 1 ,...,£ a )eCr\Cf 4=1 
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The proof of Proposition l2.4l is involved, and we postpone it to Section I2T31 
For now, we assume that Proposition 12.41 and hence Theorem 12.31 hold, and 
explain how Theorems 11.11 and 11.21 follow from this. 

Proof of Theorem By splitting the sum over the set of loop configura- 
tions C r in Theorem 12.31 according to the lengths of the individual loops, 
using (II. 6p we can write 

oo oo 1 S / \ 

s(*) = i+££i £ n( £«(**)) 

r=l s=l riH \-r s =r i=l^i£C ri 

oo oo 1 s 

M = 1+ ££^ £ II 

r=l s=l riH \-r s =r i=l 

Now suppose that, given x, there exist 7 E (0, 1) and C < 00 such that 

(2.4) \f r (x)\<CY for all r. 

For future reference, we note that this condition is implied by the stronger 
condition that 

(2.5) £ \w(£; x)\ < C 7 r for all r. 

ieCr 

Under condition (|2.4[) . if we write h(r, s) for the summand in (|2.3|) . we have 



\h(r,s)\ 



riH hr s =ri=l 7 



and thus 

00 00 00 „„ 00 



EE«^)i < E ^E (1 1 IV = ex P - 1. 

s=l r=l s=l ' r=s v 7 \ ' / 

Hence, we can apply Fubini's theorem to interchange the order of summation 
over r and s in (j2.3|) . which yields 

00 1 00 s 

^) = ! + £i£ £ n/*(*)- 



S! 

s=l r=lriH |-fs=fj=l 



Note that under condition (|2.4p . / r (x) is absolutely convergent. We now 
apply Mertens' theorem, which says that if a series Yl r ar converges abso- 
lutely, and the series Y^ r b r converges, then their Cauchy product converges 
t° (Yl r a r)(Ylr br)- I* 1 particular, by induction, the s-fold Cauchy product of 

the series ^ r a r with itself, which is ^ r 5ZrH hr s =r a riGV 2 ' ' ' a r s) converges 

to (X^r a ^') S - Applying this to a r = f r (x), we obtain 

OO^^OO \ s ✓ 00 

(2.6) Z(*) = l + £± (£/r(x)J =exp(£/ r ( i r) 

s=l S " V=l 7 V=l 

Observe that this result holds already under the weaker of the two condi- 
tions (|2.4p and f|2.5j) . but that under the stronger condition (|2.5p . the loop 
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weights can in fact be summed in any order. We will now show that the 
condition of Theorem 11.11 implies (12.5f) . Indeed, under the condition of The- 
orem [T7T|, there exists 7 G (0,1) such that (d uv — l)\x uv \ < 7 for all edges 
uv S E. Observing that if a loop takes a step along uv, then there are at 
most d uv — 1 possibilities for the next step, this implies that the sum of 
over all loops I of n steps is bounded by |F|7 n . Since a loop of 
length r takes at least r steps, summing over n > r yields (|2.5p . □ 

Proof of Theorem \1.2l Recall definition (jl.7p of the entries of A(x), which 
are indexed by the directed representative edges of G. We can interpret this 
matrix as a transition matrix for non-backtracking paths on the graph G' 
which is represented by G. This represented graph G' can be obtained 
from G by removing every chain of additional edges from G, and identifying 
the two vertices at the ends of this chain (see Figure [2] for an example) . 

Indeed, consider two directed representative edges uv and wz ^ in G, 
and write uv' and wz' for the corresponding directed edges in the represented 
graph G' . By construction, a non-backtracking path in G' can make a step 
from uv' to wt' if and only if v is linked to w in the graph G, since only then 
will v be identified with w in G' . This step corresponds to either a direct 
step from uv to wz in G (if v = w), or to a sequence of steps along the chain 
linking v to w. In either case, the matrix entry ^(x) picks up all edge 
weights and turning angles associated with these steps in G. 

We can now interpret this entry as describing the weight picked up by a 
non-backtracking walk in G' when it steps from uv' to wz 1 . Viewed in this 
way, the entry of the matrix A r (x) indexed by uv and wz is equal to the 
sum of the weights of all non-backtracking paths in G' of r steps starting 
from uv' and ending on wz'. In particular, the sum of the diagonal entries 
of A' r (x) is equal to the sum of the weights of all non-backtracking paths 
in G' of r steps starting and ending on the same directed edge. 

Now consider a loop £ of length r in G. Note that it is possible to start 
traversing this loop from each step it takes along a representative edge in two 
directions. Mapping the paths thus obtained to the represented graph G' 
yields precisely 2r/m(£) different non-backtracking paths of r steps in G' 
that start and end on the same directed edge. By (|1.2p . f)1.5|) and (|1.6p . it 
now follows that 

trA r (x) = -2rf r {x), 

where the minus sign comes from the minus sign in the definition fjl .5[) of 
the sign of a loop in terms of its winding angle. Expressed in the eigenvalues 
Xi(x) of A(x), we therefore have that 

(2-7) /r(*) = ~J>'(s). 

i 

In particular, condition (|2,4p is satisfied if maxj|Aj(x)| < 1, so in this case 
the same argument as in the proof of Theorem 11.11 yields (|2.6p . Moreover, 
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if maxj| \i(x)\ < 1, then using (12 .7p we can write 

\ r=1 i / \ j r=1 / 

and since YlT=i uV / r = ~~ in (l ~~ u ) if M < 1) we conclude that 

Z(x) = JJ(l - Xi(x)) l/2 = y/det{I - A(x)). □ 

i 

2.3. Cancellation of non-edge-disjoint loop configurations. We now 

turn to the missing step in the proofs of Theorems 11.11 and 11.21 which is 
the proof of Proposition 12.41 That is, we must show that the weights of 
all loop configurations (£i,...,£ s ) which are not edge-disjoint and satisfy 
r(£±) + • • • + r{£ s ) = r for a given r, cancel each other. What complicates 
matters here, is the fact that these loop configurations do not cancel each 
other one by one, see for example Figure l5fl Our strategy of the proof is to 
map loop configurations to so-called labelled loop configurations, which do 
cancel each other one by one, and show that this implies cancellation of the 
unlabeled loop configurations for combinatorial reasons. 

We will therefore start by introducing the notion of a labelled loop, and 
work our way from there towards the notion of a labelled loop configuration, 
and the proof of their cancellation. In words, a labelled loop is a loop with 
a label attached to each step it takes, where the labels are distinct positive 
integers. For periodic loops, the first step is repeated after completing a 
period, and we require that the label of the first step of the loop is smaller 
than the label associated with each of these repetitions. 

Formally, a labelled loop £9 is a sequence (uo, do, i>i, oi, . . . ,Un-i, ©n-i) 
satisfying the following conditions: 

LI £ = (yo, . . . , v n -i) is a loop; 

L2 {flQ, a±, . . . , a n -i) is a sequence of distinct positive integers, called 

the labelling of the loop; 
L3 if £ is periodic, i.e. m(£) > 1, then ao is smaller than af. n / m m for all 

k£{l,2,...,m(£)-l}. 

We call the number the label on step i + 1 of the loop £; we also regard 
it as a label assigned to the edge fjfj+i. We will use the superscript for 
labelled loops, and the unlabelled loop corresponding to a labelled loop will 
consistently be denoted by dropping this superscript: if £^ is a labelled loop, 
then £ is the corresponding unlabelled loop, and so on. 

Observe that one of the effects of labelling loops is that it breaks the 
periodicity of periodic loops: sequences representing labelled loops cannot 
be periodic. Therefore, if £ is periodic, we do not assign to the labelled loop 



A picture of the same configurations appears in j6 to point out the error in 
Vdovichenko's paper; it is crucial here to take the multiplicity of the loops into account. 
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1! X 3 



X 1 



X 1 X 



X 1 X 1 X 1 



Figure 5. Four loop configurations on the same vertices 
and edges, where the traversals of the same edge have been 
drawn slightly apart to make them discernible. The factors 
h lli=i sgn (-^i ) /m i ) are spelled out below each loop config- 
uration to show that the sum of their signed weights is 0. 

(vo,ao, . . . ,v n -i,a n —i) the same weight as to its unlabelled counter- 
part. Instead, we define the weights of labelled loops in general by 

ra-l 

(2.8) w(lO;x) = sgn(£)l[x ViVi+ ,, 

i=0 

where the sign is defined in terms of the winding angle of i by (|1.5|) . as before. 
Note that this weight is actually independent of the particular labelling of 
the loop, and that w(£^;x) = m(£)w(£;x). 

We write n(£) for the number of steps of a loop £ (recall that this is 
not necessarily the same as the length r{£) of the loop). By a labelled loop 
configuration we mean a collection {£f , . . . ,£f} of labelled loops, in which 
all labels are distinct and take values from the set {1,2, ... , X^i=i n (^«)}- I n 
particular, any loop configuration (£\, . . . ,£ s ) can be turned into a labelled 
loop configuration by attaching a label to every step of every loop in such 
a way, that condition IL3I above is fulfilled for every labelled loop obtained, 
and all labels 1,2,..., ^4=1 n (^i) are used. 

Now fix r and n, and consider a loop configuration (£\, . . . ,£ s ) which is 
not edge-disjoint and satisfies Ya=i r (^«) = r an d Si=i n (^«) = n - Let t 
denote the number of distinct loops in (£±, . . . ,£ s ), and write k\, . . . ,kt for 
the respective number of times each of them occurs, so that k\ + - ■ ■ + kt = s. 
Consider the collection of all labelled loop configurations {£f , ... ,£f} that 
can be obtained from {£\, . . . ,£ s ) by labelling the loops, as described above. 
For a periodic loop £t, only one of the rotations of its labelling, rotated over 
a multiple of the smallest period, satisfies condition IL3I Furthermore, in- 
terchanging the labellings of two identical loops £{ and £j yields the same 
labelled loop configuration. Therefore, the number of labelled loop configu- 
rations we obtain from (£1, . . . ,£ s ) is precisely 



We assign to each of these labelled loop configurations the same weight 
Y\t=i w {^fi x )) where we use the fact that according to the definition (|2.8p . 
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w(if;x) does not depend on the actual labelling. Then the total weight of 
all labelled loop configurations associated with (£i,.- - ,£ s ) is 

, s . s 



We claim that this is exactly n! times the total weight that all the per- 
mutations of the loop configuration (£i,...,£ s ) contribute to the sum in 
Proposition 12.41 Indeed, there are precisely 

s! 



such permutations, and the weight each of them contributes to the sum is 



1 s 

— Y[w(£i;x). 



We conclude that to prove Proposition 12.41 it suffices to show that for 
given n and r, the weights of all labelled loop configurations {if , . . . ,£f} 
such that Ylt=i n (^i) = n i J2t=i r (^i) = r an d {£■!■,■ ■■ As) is not edge- 
disjoint, sum to 0. Write Cf lT for this collection of labelled loop configu- 
rations. We will now prove the desired cancellation of weights, and hence 
Proposition 12.41 by finding a bijection g: C% r — > C%, r which maps each la- 
belled loop configuration to a labelled loop configuration which has a weight 
of the opposite sign, but with the same absolute value. 



Proof of Proposition \2.4\ Before we go into the formal details of the bijec- 
tion, let us give an informal description of how it will work. Consider a 
labelled loop configuration {£f,...,£f} G C% r , and let be the set of 
edges in G that are assigned more than 1 label in this configuration. Find 
the smallest of all the labels that are assigned to the edges in E^ , let a be 
this label, and let uv be the edge to which this label is assigned. Next, find 
the second smallest label b which is assigned to the edge uv. 

The label a labels a step of one of the loops ii. The label b either labels 
another step of the same loop £i, or it labels a step of a second loop £j, i ^ j. 
The bijection involves interchanging the "connections" on one side of the two 
steps marked a and b (either at the vertex u or at the vertex v ), as illustrated 
in Figured! It is clear that this operation does not change the absolute value 
of the weight of the configuration, since the total number of steps that go 
through a given edge does not change. But Figure [6] also suggests that the 
operation corresponds to increasing or decreasing the number of "crossings" 
in the configuration by 1, which should indeed lead to a change in sign. 

However, signs were formally defined in terms of winding angles, not 
numbers of crossings, since it is more difficult to make sense of the latter 
when loops are not edge-disjoint. Furthermore, we must still formally define 
the mapping g. We will now deal with these technical issues. 
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Figure 6 . All cases that occur in the cancellation of labelled 
loop diagrams, as explained in the text. The curves pi and p2 
represent arbitrary paths connected to the vertices u and v. 

For the formal treatment of the bijection, we need to introduce some 
additional notation. Given a sequence a = (ao, . . . , a n ) of arbitrary elements, 
we write a~ l for its reversion a -1 = (a„, a n -\, . . . , ao). lib = (bo, • • • , b m ) is 
another sequence of arbitrary elements, we write a© b for the concatenation 
of a with b, that is, 

a© b = (ao, ...,a n ,b ,... ,b m ). 

The weight of a labelled loop configuration {if , ■ ■ ■ ,£f} is defined as the 
product of the signs of the loops l\, . . . ,£ s , times the product of all the 
edge weights picked up by all the loops. As was anticipated above, the 
product of edge weights will not change under the bijection, so we will only 
be concerned with the product of the signs of the loops. We recall from (II. 4p 
and (II. 5p that the sign of a loop £ = (vq, . . . , u n -i) is defined in terms of its 
winding angle as 

(2.9) sgn(t) = -exp(~a(£)), 

where the winding angle a(i) is given by 

n-l 

(2.10) a(£) = ^2 Z(v i+1 - Vi,v i+2 - v i+1 ). 

i=0 

We now define the winding angle and sign of a rooted loop (vq, . . . , u n -l) 
by the exact same formulas. In particular, all rotations of a loop £ have the 
same winding angle and sign. On the other hand, the reversion t^ 1 of I and 
all its rotations are traversed in the opposite direction, and therefore they 
all have winding angle a(£ _1 ) = —a(£). However, since the winding angle 
of a loop is a multiple of 2tt, we do have that 

(2.11) sgn^" 1 ) = sgn(£) for all rooted loops £. 
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We call all the rotations of a loop i, and all rotations of its reversion l~ , 
alternative representations of I. All these representations have the same 
sign. Likewise, the rotations of a labelled loop i^ and its reversion v^~ x will 
be called representations of this labelled loop. 

We also need to define the winding angle for paths in G which are not 
loops. Recall that a path p = (vq, . . . ,w n -i) is not a loop if UoU n -l ^ E, 
vo = v n -2, or v\ = v n -\. If we follow such a path from vq to we turn 

through n — 2 angles, and it is natural to define the winding angle of p by 

n— 3 

We now have all the notation we need to define and analyse the bijection 
formally. So consider a labelled loop configuration {if , . . . ,if} 6 C% r , and 
define , a, b and as above. We will now explain to which labelled loop 
configuration our configuration {if , . . . ,if} is mapped by the bijection, and 
prove that the image has the opposite sign, and hence the opposite weight. 
There are three possible cases to consider, which are illustrated in Figure El 

Case 1: The labels a and b belong to different labelled loops. Let if 
be the labelled loop containing label a, and let i'f be the labelled loop 
containing label b. Then these labelled loops have representations of the 
form if = (u,a,v) © pf and £J = (u,b,v) © pf , respectively, where pf 

and pf are paths interspersed with labels. We can now form the combined 
representation 

Ifj = (u, a, v) © pf © 0, b, v) © pf 

of a new labelled loop ifj. Our bijection maps {if, . . . ,if} to the labelled 
loop configuration 

{£f,...,if,£%}\{£?,£f}. 

To see that this labelled loop configuration has the opposite sign of its pre- 
image {if , . . . , if}, note that by (p^fpTTl) . 

sgn(^) sgn(^) = sgn(ij) sgn(lj) = expQa(ij) + ^.{Ifyj 

= expQa(ly)^ = -sgn(ly) = -sgn(4,-). 

Case 2: The labels a and b are on steps of the same labelled loop taken 
in the same direction. This case is the reverse of Case 1. The labels a and b 
are in a labelled loop if which has a representation of the form 

if = (u, a, v)®pf® (u, b, v) © pf . 

From this we obtain the representations (u, a, v) © pf and {u, b, v) © pf of 
two new labelled loops if and if. The bijection maps {if , . . . ,if} to the 
labelled loop configuration 

{if,...,if,if v £f 2 }\{if}. 
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The same argument as in Case 1 shows that sgn(^i) sgn(^2) = — sgn(^). 

Case 3: The labels a and b are on steps of the same labelled loop taken in 
opposite directions. In this case the labels a and b are in a labelled loop if 
which has a representation of the form 

if = (u, a, v) © pf © (v, b, u) © pf. 

From this we can construct the representation 

P = o, a, v) © pf - 1 e (v, b, u) © pf 

of a new labelled loop I®. The bijection maps {if, . . . ,if} to the labelled 
loop configuration 

{if,...,if,iO}\{if}. 
To verify that these loop configurations have opposite signs, observe that 

(2.12) a(ii) = a((u, v) © p x © (v, u)) + a({v, u) © p 2 © (u, v)) , 
and likewise 

(2.13) a(t) = a((u, v) © pj" 1 © {v, u)) + a({v, u) © p 2 © (u, v)) , 

where pi and p2 are the paths obtained from p^ and pf by dropping the 
labels. Now notice that upon reversion, 

(2.14) a((u, v) © pi © (v, u)) = -a((u, v) © pj" 1 © (v, u)) . 
Furthermore, it is not difficult to see that 

a((u, v) © pi © (u, n)) = 2m7r + it for some m G Z. 
Together with (I2TT21) . (l2~T3l) and pii|> . this implies 

sgn(£j) sgn(^) /« ~ i 



expQa^) - = -1. 



sgn(£) S gn(^) 

We conclude that in all cases, the labelled loop configuration {if , . . . , if} 
is mapped to a labelled loop configuration of opposite weight. From the 
explicit descriptions given above, it is not difficult to see that the mapping 
is bijective. As we have explained above, this implies Proposition 12.41 □ 

3. Proofs of our results for the Ising model 

In this section, we will apply Theorem 11.21 to the Ising model on the 
square lattice 1? . We will see that this leads to explicit expressions for the 
free energy density and two-point functions in terms of sums over loops in Z 2 
or its dual Z 2 *, valid all the way up to the critical point. To get started, we 
must find a bound on the spectral radius of the matrix A(x) appearing in 
Theorem 11.21 This will be done in Section 13.11 Then we will study the free 
energy density in Section T3.21 and finally the two-point functions at low and 
high temperatures in Sections 13.31 and 13. A\ respectively. 
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3.1. Bound on the spectral radius. Let G = (V,E) be a fixed finite 
rectangle in with no additional edges (i.e. Ea = 0)- Without loss of 
generality, we may assume that the vertex set is 

(3.1) V = {0,1,...,M-1} x {0,1,..., N -1}. 

Since we are on the square lattice, directed edges can point in only 4 direc- 
tions, and we now introduce some convenient notation for this specific case. 
We write vl, vl, v-* and u«- for the directed edges from v to, respectively, 
v + (0, 1), v— (0, 1), v + (l, 0) and v — (1, 0). We also write \v for the directed 
edge pointing from v — (0, 1) to v; iv, -*v, <-u are defined analogously. 

Given a vector of edge weights x = (x uv ) uvt =E on E, A(x) is the transition 
matrix indexed by the directed edges of G, as defined in (jl.7p . For a vertex v 
not on the boundary of G, the row of A(x) indexed by for instance, has 
exactly 3 nonzero entries, corresponding to the 3 possible steps that a loop 
can take from -*v. To be precise, with u = v — (1,0), these 3 entries are 

A- iv^v-* (x) — %uvi -A— >v,vt (p^) — %uv& ^ ; *v,v\, (p^) — %uv& ^ • 

Observe that most rows of A(x) have exactly 3 nonzero entries. The only 
exceptions are the rows indexed by directed edges pointing to a vertex 
in dG. These exceptional rows make it impossible to compute the eigen- 
values of A(x) directly. We will therefore wrap G around a torus, so that 
all vertices can be treated alike, and then bound the eigenvalues of A(x) 
in terms of those of the transition matrix on the torus, or equivalently, the 
transition matrix for a graph in the plane representing the torus. 

To be precise, we first extend our graph G to a graph G® (© stands for 
"torus"), by adding edges and vertices as shown in Figure [2] (left). Note that 
this adds directed representative edges v-* and *-v for every vertex v on the 
right boundary of G, and v\ and lv for v on the top boundary. All other 
edges that are added are considered as additional edges in the graph G®. 
Henceforth, when we work on the graph G®, computations will be performed 
modulo M and N in the two respective lattice directions. 

The matrix A® is the transition matrix for the graph G® , with specific 
edge weights chosen as follows: all representative edges of G® have edge 
weight 1; for the additional edges, we choose the edge weights in such a way, 
that the product of the edge weights along every chain of additional edges 
linking opposite sides of the rectangle to each other, is —1. Note that by 
this choice, the factor —1 will exactly compensate the sign picked up by a 
path which follows the chain, because of the 4 quarter-turns it makes. 

Proof of Theorem \1.3X We first prove that the spectral radius of the ma- 
trix A® is v2 + 1. To this end, assume that the rows of A® are arranged in 
such a way, that for every vertex v G V , the 4 rows indexed by -*v, lv, <-u 
and iv immediately succeed each other in this order. Let II be the permu- 
tation matrix which permutes the columns of A® so that column vd maps 
to column dv, for all v E V and d 6 {t, i, <-}. 
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By construction, the matrix A®LT with the permuted columns is now 
a block-diagonal matrix, since the 4 rows indexed by the directed edges 
pointing to v are matched along the diagonal with the 4 columns indexed by 
the directed edges pointing out from v. By considering the turning angles, 
it is easy to see that each 4x4 block is equal to the Hermitian matrix 

1 exp(i7r/4) exp(— i7r/4) 

exp(— in/4:) 1 exp(i7r/4) 

exp(— z7r/4) 1 exp(i7r/4) 

exp(i7r/4) exp(— wr/4) 1 

which has eigenvalues y/2 + 1 and \/2 — 1, both of multiplicity 2. 

For a square matrix B of dimensions d, let \\B\\ denote the operator norm 
induced by the Euclidean norm on C . Since A is Hermitian, its spectral 
radius is equal to \\A\\. It follows that the operator norm of A® IT is given by 
||A|| = y/2 + 1, and since permuting columns does not change the operator 
norm of a matrix, we conclude that ||A®|| = y/2 + 1. 

In general, the spectral radius of a matrix B is bounded from above by 
the operator norm [|J3||. We will now use this fact, together with the sub- 
multiplicativity of the operator norm, to bound the spectral radius of the 
matrix A(x). To this end, let D(x) be the diagonal matrix of the same 
dimensions as A®, defined as follows. For vertices v on the right boundary 
of G, the diagonal entries of D(x) on the rows v-> and *-v are 0, and so are 
the diagonal entries on the rows and iv for v on the top boundary of G. 
For all other directed edges uv in the graph G®, the diagonal entry of D(x) 
on row ui is equal to the edge weight x uv . 

Now consider the matrix D(x)A®. This is just the matrix A® with all 
rows corresponding to directed edges uv in the graph G multiplied by x uv , 
and all rows corresponding to directed edges which are in the graph G®, 
but not in the graph G, zeroed out. Therefore, considered as a transition 
matrix on the graph G®, D(x)A® generates exactly the same loops as the 
matrix A(x), with the exact same weights. 

Let ij) be an eigenvector of A(x), indexed by the directed edges in G. 
Extend this vector to a vector ip® indexed by the directed representative 
edges in G® , by setting ijj®^ = ip^ if uv is an edge in G and ip®^ = if uv is 
a representative edge in G® but not in G. Then, by construction, tp® is an 
eigenvector of D(x)A® with the same eigenvalue as ip. Hence all eigenvalues 
of A(x) are eigenvalues of D(x)A®, and we conclude that 

max|A;0)| < ||D(x)A®[| < \\D(x)\\ ■ \\A®\\ = (V2 + ljlMI^. 

i 

The desired bound on |/r(^)| now follows from (|2.7p and the fact that the 
number of directed edges in G is bounded by 4|V|. □ 

Remark. Using Fourier transforms, we can actually compute all eigenvalues 
of A®. Also, with some extra effort, it is possible to show that for all finite 
rectangles, the spectral radius of A(x) is strictly less than (v2 + 1) ll^llco- 
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3.2. Free energy density. We are now going to use the bound obtained 
in Theorem 1 1.3 1 to derive an explicit expression for f((3) as a power series in 
the parameter exp(— 2/3) for /3 > /3 C , and in the parameter tanh/3 for /3 < /3 C . 
From this expression we then also obtain Onsager's formula. 

The coefficients of the power series for are determined by the loops 
in 1? for which the origin o = (0,0) is the lexicographically smallest vertex 
they traverse. To be more specific, we define C° as the collection of all loops 
i = (vo, . . . , v r -i) in 1? such that vq = o. Here, all edges in 1? are regarded 
as representative edges, so that the length r of a loop is equal to its number 
of steps. Let all edges of 1? have the same edge weight x. Then we can 
define the weights w(£;x) of all loops I G C° by (jl.2p . as before (by slight 
abuse of notation, x now stands both for the weight of a single edge, and 
for the vector of all edge weights). Write 

Observe that f°(x) = f°(l)x r , so ^ r f°(x) is really a power series in x. 
Theorem 3.1. For all (5 G (0, /3 C ) U (/3 C , oo), the free energy density satisfies 

oo 

(3.2) - /3f(P) - ln(2 cosh 2/3) = - ln(l + x 2 ) + ^ f°{x), 

r=l 

where x = exp(— 2/3) for j3 G (/3 C , oo) and x = tanh /3 for /3 G (0, /3 C ). 

Proof. We start with the high-temperature case, so fix /3 G (0, /3 C ) and set 
x = tanh /3. Let G be a rectangle in Z 2 , and take the set of additional edges 
to be empty. Note that by (Oil) , x G (0, \/2 — 1). By (flTTSl) . we have 

(3.3) = |^|ln2 + |-E|ln(cosh/3) +lnZ G (a;), 

where Zq{x) is the generating function for the graph G with edge weights 
equal to x. By Theorems 11.21 and 11.31 In Zq{x) is equal to J2 r fG,r(x), where 
fc,r{x) denotes the sum of the weights of all loops of length r in the graph G. 

Consider these loops of length r in G. For each vertex v G V, let C^{G) 
denote the collection of loops in G of length r for which v is the smallest 
vertex traversed. Observe that if v has distance at least r to the boundary 
of G, then C^.(G) can be mapped bijectively to C° by a translation on Z 2 , 
hence J2ieC v (G) w ^'-> x ) = f°( x )- There are at most \dG\r vertices at a 
distance less than r from dG, and since \x\ < 1, for such a vertex v we have 
YleeC v (G)\ w (^> x )\ — 3 r (by counting non-backtracking paths). From these 
observations and the fact that |<9G|/|y| — > as G — > Z 2 , it follows that 

lim -— fG,r( x ) = lim 2 ]T7\Z^j w (^ x ) = fr( x ) for all r > 1. 

11 11 vevt£Cy(G) 
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Furthermore, Theorem 11.31 says that |/c, r (af)| < 2\V\{V2 + l) r x r . Therefore, 
by dominated convergence and Theorem 11.21 



oo 



W\ Zg{X) = ^ W\ hAx) = ^ ^ 

r=l r=l 

We now combine this with (|3.3p . and use lim G _ >Z 2|.E7|/|V'| = 2 to obtain 

oo 

-/3/(/3) = G hm 2 — In = ln(2 cosh 2 (3) + £ f°(x). 

Using cosh 2 /3 + sinh 2 /3 = cosh 2 /3(1 + x 2 ) = cosh 2/3, we can cast this into 
the form of equation f|3.2|) . We perform this final step here to obtain a 
common expression for /(/?) for all non-critical temperatures. Indeed, the 
low-temperature case can be treated in a similar manner, except that one 
must work on the dual graphs G* with edge weights x = exp(— 2/3) on the 
dual edges, and use (|1 . 12[) instead of fll . 13 j) . In this case one obtains 



1 oo 
r=l 



which can also be cast into the form of equation (|3.2p . by using the logarithm 
of the equality 2 cosh 2/3 = e 2/3 (l + x 2 ) to rewrite the 2/3 term. □ 



Proof of Theorem \1.4\ The expression for /(/3) derived above establishes di- 
rectly that the free energy density is analytic on (0, /3 C ) U (/3 c ,oo). We will 
now show that /(/3) is given by Onsager's formula. 

In the proof of Theorem I3.ll we have obtained ^2 r f°(x) as the limit of 
l^l -1 fG,r{x)- It is clear from the proof that here we may as well replace 
fc,r(x) by the corresponding sum of loop weights for the periodic graph G® 
from Section \3. 11 In fact, the argument becomes even simpler on G®, since 
we no longer have to treat vertices near the boundary separately. The tran- 
sition matrix generating the loops in G® with the desired edge weights x 
is xA®. Hence, by Theorem 11.21 and the proof of Theorem II .31 in Section 
we see that for x G (0, \Jl — 1), 

oo 1 oo 1 1 

(3-4) £ = Jim — £ fa. tr (x) = Jim — - lndet(I - xA®). 

r=l r=l 

We can compute det(I — xA®) by taking the Fourier transform of A®, a 
computation which has appeared before, for instance in [9ll21}[25]. We give 
the calculation here for completeness, but will go through it a bit quickly. 

Without loss of generality, we may assume that V is the set (|3.ip , in which 
case the Fourier transform of A® is defined as 

M-l JV-1 

X® - _L_ v V P -¥(p k ~p' k ')-¥(i l -i' 1 ') a® 
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where d, d! G {t, I, <-}• The calculation of this Fourier transform is made 
straightforward by the periodicity of A®. The calculation reveals that the 
only entries surviving the summations are those for which p' = p and q' = q. 
Hence, A® is a block-diagonal matrix of 4 x 4 blocks. To be precise, writing 
u p = 2irp/M, u q = 2-irq/N, the 4x4 block for given p and q is 



A® 



J.uip+iir/A q giLUp-iTT/4 

q gio; ,2+171-/4 Q 



e iuj q —in/4 e iui q 

q e —iu p — in/4 e ~iui P e —iui p +in/4 

p— iu q +in/4 q ^—iujq—in/4 ^—iu q 



Since det(I — xA®) = det(I — xA®), from this Fourier transform we obtain 
M-l N-l 



(p,q)-,(p,q) 



det(I -xA®)= Yl "Q det V ~ xA l 

p=0 q=0 
RI-l N-l 

= ] [ Y\ [(1 + x 2 ) 2 - 2x(l - x 2 )(cosWp + cosa; g )] . 

p=0 g=0 



Using (|3.4p . we conclude that 

oo _. 

(3.5) V/°(x)= lim — — lndet(I-xA®) 



r=l 



1 



2n r-2n 



„2\2 



2x(l — x 2 )(cosu;i + COSW2)] du\ du2- 



To finish the computation, note that by (|3.2p . we have 
~2 cosh 2/3" 



-/3/(/3) = In 



1 + x 2 



r=l 



1 

8^ 



/ / ln 


"4 cosh 2 2/3" 


/0 •/ 


(l + x 2 ) 2 _ 



r=l 



Combining this with (|3.5p . and then using the identity 

2x(l - x 2 ) _ sinh2/3 
(1 + x 2 ) 2 ~ cosh 2 2/3' 

which holds both for x = exp(— 2(3) and for x = tanh/3, we obtain 

— /5/(/3) = - — 2 / / ln[4cosh 2 2/3 - 4sinh2/3(cos w x + cosu 2 )] du\ du 2 . 
8 71 " Jo Jo 



This is Onsager's formula for the isotropic Ising model on 



□ 
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3.3. Low-temperature correlations. In this section, we discuss the Ising 
model with positive boundary conditions. We consider rectangles G = 
(V, E) in 1? (which later tend to 1?) and denote by G* the weak dual 
of G. Recall that every spin configuration a £ S7 + on G corresponds bijec- 
tively to an even subgraph of G*, that is, a graph in which all vertices in V* 
have even degree. For given a, we denote the corresponding even subset 
of E* by F{a)\ for given even F C £*, we denote the corresponding spin 
configuration by o~(F). 

Setting x e = e~ 2 @ for every edge e in 1?* , by (jl.lip and ()1.12p we have 

(3.6) P ^ {(T) = -}— yi Xe , aen +, 

G V ' e€F(a) 

where Zq* (x) is the generating function for G* with edge weight vector 
x = {x e ) e& E*- Note that here we have defined the edge weights on all edges 
of Z 2 *, which is necessary because our graphs G will tend to Z 2 . However, 
the weight vector x is implicitly restricted to the relevant edges when we 
work on a given graph G. Such implicit restrictions to the relevant edges 
will occur throughout this and the following section. 

Proof of Theorem \1.5\ Fix u,v G 1? , u ^ v, and let 7 be a self-avoiding 
path in 1? from u to v. We may assume that G is large enough so that u, 
v and 7 are all contained in the area spanned by G, see Figure [3 We will 
express the two-point function (a u a v )^, „ as the quotient of two generating 

functions. To this end, we define new edge weights x' e on the edges of 1?* 
such that x' e = —x e if e crosses 7, and x' e = x e otherwise. The reason for 
defining the weights x' e in this way is the crucial fact that for all a G 

(3.7) G U G V Yl x e = Yl X 'e- 

e£F(a) eG-F(cr) 

To see this, recall that the edges in F(o~), by their very definition, cross 
edges xy £ E for which a x / a y . If a u = a v , then following 7 from u to v , 
we necessarily cross an even number of such edges. If a u 7^ o~ v , then we cross 
an odd number. In either case, (13.7P holds. 
With the help of (|3.6h . we can write 

(VuO- v )g,i3 = X] VuVvPg^o-) = z \s Yl GuCJv II Xe ' 
aen+ G *^ x ' uen+ eeF(a) 

and using (|3.7|) and the bijection between even F and VL + , we obtain 

(3 . 8) Mh=^ } e n 

^ v ' even FcE* egF ^ v ; 



The idea that correlations in the Ising model can be studied by means of 
ratios of generating functions with changed edge weights (or equivalently, 
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Figure 7. The paths 7 (with bold edges), that we use to 
study the 2-point functions {0 u ct v )q p. The low-temperature 
case is on the left (spins are —1 in the gray squares, +1 in 
the white regions), the high-temperature case on the right. 

changed spin-spin interactions) has arisen before in the physics literature, 
see |13| . It now follows from Theorems 11.21 and 1 1 . 3 1 that for f3 > f3 c , we have 

(cr u o v )^p = expf ^2 ^2 [ w (^ x ') ~ w (ti x )]h 
V=i eec r (G*) ' 

where C r (G*) is the collection of loops of length r in the graph G*. 

Call a loop in G* uv-odd if it crosses 7 an odd number of times. It is not 
difficult to see that this notion depends only on u and v, and not on the 
particular choice of 7. Observe that for uv-odd loops £, w(£; x') = —w(£; x), 
while for loops £ that are not uv-odd, w(£;x') = w(£;x). It follows that 




l MD-odd 



Note that a uv-odd loop of length r cannot travel far from u and v. To 
be precise, these loops must be contained in U B%, where B™ is a square 
in the plane of side length r centred at u, and B^ is defined similarly. To 
study the convergence of (|3,9p as G — > 1?, for arbitrary rectangles R in 1? 
that can be finite or infinite, and even equal to 1?, we now define 

a r (R*;x):= w(£;x). 

l uv-odd 

This definition makes sense both for finite and infinite R, since the loops 
contributing to the sum must be contained in B™ U B^.. 

Let B™ denote the smallest rectangle in R 2 containing both B™ and B%, 
and write R*nB^ v for the largest subgraph of R* which is a rectangle in T?* 
entirely contained in B^ v . Then for all R, 

(3.10) a r (R*;x) = a r (R*nB™;x) = - ^ [w(£; x) - w(£; x')] . 
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Now, since the volume of B™ v is bounded from above by (\\u — v\\ +r) 2 , and 
exp(2/3 c ) = y/2 + 1 by (|1.14p . Theorem 11.31 yields the uniform bound 

(3.11) \a r (R*; x)\ < 2(||u-u|| + r) 2 exp(-2(/3 - /3 c )r) for all R. 

We now return to §3M). Since eventually, G* n B? v = Z 2 * n B?° when 
G — > Z 2 , from (|3.10p we conclude that 

a r (G*;x) — > a r (Z 2 *;x) for all r > 1. 

Moreover, the a r (G*;x) are uniformly bounded in G by the right-hand side 
of (|3.1ip . which is summable over r. Therefore, by dominated convergence, 

oo oo 

lim y~^a r (G*;x) = a r (Z 2 *; x) , 

^ r=l r=l 

where the series on the right is absolutely summable. Using (|3.9p . this gives 
the following result, which we state here as an intermediate theorem because 
of its significance, and for future reference: 

Theorem 3.2. For all (3 G (fi c , oo) and fixed u, v G Z 2 (u ^ v), 

Yaa n {a u o- v )% g = expl -2 ^ ^ io(^;a;) ) =: (<r u (r„)+ 2 * > 0. 



= l£G£ r (Z 2 *): 
ut>-odd 



Next, we consider {o~ u )~q a for n G G\dG. We can treat this like {o- u o~ v )~q a 
by taking w on the boundary of G, since then = +1. We now call a loop 
which crosses 7 an odd number of times u-odd, since this notion depends 
only on u, not on v. The box B^ v can be replaced by f?" in the argument, 
which replaces (\\u — v\\ + r) 2 by r 2 in ()3.1ip . In this case, the result is: 

Theorem 3.3. For all f3 G (/3 C , 00) and /ixed u G Z 2 ; 

Jim 2 ((r u )+^ = expf-2^ 2 w(4x)J =: M+ 2)/3 > 0. 

^ r=l leC r (Z 2 *): 
i u-odd 

The limit in Theorem 13.31 is easily seen to be independent of the choice 
of u, and we take u = o as the canonical choice. We now consider what 
happens to the two-point function when we take u and v further and further 
apart. When r < \\u — v\\/2, the boxes B™ and B% are disjoint. If this is 
the case, a uv-odd loop of length r in 1?* must be contained in either B™ 
or B^.. Hence we have 

(a u a v )^ 2 ^ = expf -2 ^ a r (Z 2 *;x)\ 

^ r>||«-i;||/2 

xexp (-2 [a r (Z 2 * nB?;x) + a r (l?* nB*;x)] J. 

^ r<||«-t;||/2 
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When \\u — v\\ —> oo, the first factor converges to 1 exponentially fast, since 
the uniform bound in ()3.11j) applies to a r (Z 2 ;x). In the second factor, 
the first term in the sum is a sum over the u-odd loops of length r, and 
the second term is a sum over the f-odd loops. Hence the second factor 
factorizes and converges (exponentially fast) to [(0" o )i 2 J 2 . □ 



3.4. High-temperature correlations. In this section, we discuss the Ising 
model on rectangles G = (V, E) in 1? (which will again tend to 1?) with free 
boundary conditions. Our aim is to study the two-point functions {(T u (^v)g% 
at high temperatures. Prom the definitions (jl.9p . (jl.lOp and (|1 . 15[) we have 

<re£"2 tree G '@ cref2 free xyGE 

Performing the high-temperature expansion on the right-hand side of this 
expression, in the way explained in Section 11.2.11 leads to 

free 2l^l(cosh/^l ^ n 

G >P FCE: 5F={u,v} eeF 

where we set x e = tanh /3 for every edge e in 1? , and 5F denotes the set of 
all vertices that have odd degree in (V,F). Using (|1.13[) . we conclude that 

(3.12) Mw=zk* £ n**. 

K ' FCE: SF={u,v}e€F 

where Zq{x) is the graph generating function for the graph G with edge 
weight vector x = (x e ) eg £. 

Proof of Theorem \1.6l As in the low-temperature case, we first fix u, v E 1?, 
and we will again work with a path connecting u with v. However, 
this time, we will also use this path to construct a modified graph, in which 
the edges along the path are included as additional edges. 

To be precise, we first choose dual vertices u* and v* of such that 
\\u — u*\\ = \\v — v*\\ = 1/2, and then we pick a self-avoiding path 7 in Z 2 
from u* to v*. Let denote the set of vertices in 7, and let £L denote the 
set of edges traversed by 7, together with the edges uu* and vv*. For an 
arbitrary rectangle R in 1? (either finite or infinite) containing u and v, we 
denote by i? 7 the graph obtained from R by adding all vertices in to its 
vertex set, and all edges in to its edge set. In i? 7 , all edges added from 
the set Ely are considered as additional edges, and the edges from R are 
considered as representative. See also Figure [7] (right). We can in particular 
augment Z 2 in the way described above to obtain the graph Z 2 . 

As in the low-temperature case, we now define weights x' e on the edge set 
of Z 2 such that x' e = —x e if e crosses 7, and x' e = x e otherwise. We also 
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define edge weights x'(t) e on the edge set of Z 2 as follows: 

x' e if e is an edge of Z 2 ; 



^7 (*)e 



1 if e € Ey \ {uu*}; 
t if e = uu*. 



To motivate this definition, consider a given rectangle G = (V,E) in 1?, 
large enough so that u, v £ V. We claim that 

(3.13) e n^= e n 

FCE: egF even Fc£UE 7 : e£F 

5F={u,v} FdE~/ 

To see this, note that we can bijectively map every F contributing to the 
first sum to a subgraph in the second sum, by taking the union F U E~. 
Doing this may introduce edge crossings, whence the factor (—1) Cf , but 
these are compensated by switching from the edge weight vector x to x 7 (l). 

The crucial step is now to recognize the last expression as the derivative 
of a graph generating function. Indeed, a simple consideration shows that 

e (-ipiKm* = e M) Cr n <«).). 

even FcEUB 7 : egF ^even Fc£U-E 7 egF ' 

evaluated at any i, since any even F C E~ contributes at most one factor t 
to the product of edge weights on the right. In particular, we are allowed to 
evaluate the derivative at t = 0. By f)3. 12|) and (|3.13p . this establishes that 

(3.14) Z G (x) ■ {a u a v )%^ = ^Z Gi (x' y (t)) 



t=o 



We now fix p e (0, p e ), so that by (fTTTil) . x e G (0, y/2— 1) for every e£fi, 
and by Theorem II .'6\ the spectral radius of A(x') is strictly less than 1. We 
now need a similar bound on the spectral radius of the matrix A 7 (x^(t)), 
which is the transition matrix for the modified graph G 7 with edge weight 
vector x'^(t). The difference between these two matrices is that A 7 (x 7 (t)) 
allows transitions between u and v along the chain of additional edges in E^. 
This means that the 32 matrix entries from du to vd! and from dv to ud 1 , 
with d, d' 6 {t, 4, -►, <-}, are nonzero in A 7 (x' (t)j for t 7^ 0, while they are 
in A(x'); all other entries of the two matrices are the same. 

The 32 deviating matrix entries are all of the form ie*^/ 2 , where <fi is a 
sum of turning angles. Here, t will be treated as a complex variable. For 
t = 0, A 7 (x 7 (t)) = A(x'). Since the eigenvalues vary continuously with t, 
we conclude that there exists e > such that for all t satisfying \t\ < e, the 
spectral radius of A 7 (x 7 (i)) is bounded from above by some a G (0, 1). 

Hence, if |t| < e, Theorem 11.21 applies . and we obtain 



(00 
Em*) 
r=l 
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where 

f lr (t)= W&ltyW)- 

Note that, this last sum being finite, the f-y r (t) are polynomials in t. Also, 
from (|2.7p it follows that \f-y r (t)\ < 2|y|a r . Therefore, the partial sums of 
the series ^ r f^ r {i) are uniformly convergent for \t\ < e, and the sum of the 
series is an analytic function of t. Moreover, the derivatives of the partial 
sums also converge uniformly to the derivative of the sum of the series. 
From all this, it follows that the right-hand side of (|3. 14[) is equal to 



CO r, 

o 



dt ^ 

1 £e£r(G 7 ) 



t=o. 




ex p E E 




In the first factor, the only loops that survive the differentiation are those 
that visit the edge uu* , since only they contribute a factor t to the weight. 
Taking the derivative at t = 0, we are only left with those loops that visit the 
edge uu* exactly once. In the second factor, because we set t = 0, the only 
loops that contribute are those that do not visit uu* . This leaves precisely 
all loops in the graph G. The right-hand side of (|3. 14[) therefore becomes 

(E E to(/;<(l)))exp(f; Y "(V)). 

\r=l £g£ r (G 7 ): / \r=l i&C r (G) J 

I visits uu* once 

from which, applying Theorem 11.21 again to the second factor, we find that 

(3.15) K<7«)^ = 

\r=l 

I visits uu* once 

Note the ratio of graph generating functions in (13.150 . Recall that we have 
seen such a ratio of graph generating functions before in the low-temperature 
case, namely in ()3.8f) . Thus, this ratio can be interpreted as a two-point 
function between the spins at u* and v* in a dual Ising model with positive 
boundary conditions at the dual low temperature (3*, given by exp(— 2(3*) = 
tanh/3. Using (13. 9ft . we can express this ratio in terms of a sum over all 
u*v *-odd loops in the graph G, if we like. Hence, it is possible to express 
{ a u a v)<3% completely in terms of loops. 

Next, we want to consider the limit as G — > 1? . By the argument given 
in Section 13.31 we already know that the ratio of graph generating functions 
in (I3.15|) converges to (<J u *a v *)^ 2 * „„. It remains to consider what happens 
to the sum over the loops that visit uu* once. To this end, for a general 
finite or infinite rectangle R in 1? containing u and v, we define 

(3.16) a r (Ry;x' 7 ):= Y M^), 

£ visits uu* once 
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where we have simplified the notation somewhat by writing x'^ := x' (1). 

As in the low-temperature case, the loops that contribute to a r (.R 7 ;x 7 ) 
must be confined to the box B^ v , defined in the same way as before, except 
possibly for the steps taken along the additional edges in £L, which do not 
count for the length of the loop, and are allowed to go outside B^ v . Hence, 

(3.17) a r {R 1 ;x' 1 ) = a r ((RC\B™) 1 -x' 1 ), 

where R(~) B™ is the largest subgraph of R contained in B™ v , as before. It 
follows that a r (Gj)x') — > a r (Z^;x^) for all r > 1. As before, we now want 
to use dominated convergence to prove that 

oo oo 

(3.18) lim ^a r (G 7 ;x^) = ^a r (Z^;x^), 

r=l r=l 

and that the right-hand side is absolutely summable. This requires an ap- 
propriate uniform bound (in R) on the right-hand side of A3. 16f) . 

To obtain this bound, by (I3.17D it is sufficient to consider an arbitrary 
finite rectangle R in 1? containing u and v , and such that R is contained 
in B™ v ■ Note that every loop of length r in i? 7 which visits the edge uu* 
once, has a representation of the form p © 7, where p is a path of length r 
in R from v to u. Here, we use the facts that the part p of the loop never 
visits uu* , and that the steps taken along 7 from u to v do not contribute 
to the length of the loop, since they are along additional edges. 

Let Ar(x') be the transition matrix for the graph R with edge weights x' e . 
Note that the sum of the weights of all paths p of length r from v-* to tu, 
for instance, is given by the entry of the matrix A^(x') in row v-* and 
column \u. To compute the sum of the weights of the corresponding loops 
p©7, we only need to multiply this entry by the factor , where (f> is the 
sum of the turning angles encountered in the path from \u to t>-> along the 
edges in E~. From these observations, we can conclude that 

a r (i? 7 ;x 7 )<16||A^xO|| max , 

where ||-|| max denotes the maximum-entry norm, and the factor 16 comes 
from the fact that there are 4 directed (representative) edges pointing out 
from v, and 4 pointing to u. By the proof of Theorem 11.31 and the fact that 
the maximum-entry norm of a matrix is bounded by the operator norm, and 
using that (tanh/3 c ) _1 = \/2 + 1 by (|1.14p . we obtain 

(3.19) a r (R,;x' 7 ) < 16\\A r R (x')\\ < 16||A R (x')ir = 16 {^^J ■ 

We emphasize that this bound holds uniformly for all finite and infinite 
rectangles R containing u and v. Hence, (|3.18|) holds by dominated conver- 
gence, and we have established the following intermediate result: 
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Theorem 3.4. For all j3 £ (0, f3 c ) and fixed u,v G 1? (u^v), 

V=l feC r (Z2) : 

£ visits uu* once 

As an aside, we note that this result simplifies when u and v are on the 
same face of 1? (i.e. ||u — t>|| = 1), since then we can take u* = v*, so that 
o~ u *°~v* = 1- Moreover, since the path 7 is void in this case, none of the edge 
weights on G need to be modified, that is, x' = x in this special case. 

To conclude the proof of Theorem 11.61 we consider what happens to the 
two-point function when we let ||u — v\\ tend to infinity. Since the loops that 
visit uu* necessarily have length at least \\u — v\\, we can write 




By Theorem 11.51 the two-point function on the right is bounded between 
and 1. Alternatively, at this stage we could also observe that the ratio of 
graph generating functions in ()3. 15j) is always between —1 and +1, since the 
same even subgraphs contribute to both generating functions, but only in 
the numerator, some of them have a negative sign. Furthermore, the bound 
in (I3.19P holds for a r (Z^;z'). Hence we have the explicit bound 




which tends to exponentially fast as \\u — v\\ — )■ 00. □ 
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